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We study magnetic and critical properties of the alternating spin antiferromagnetic Heisenberg 
chain with S — 1/2 and 1 in a magnetic field at T = 0. The numerical diagonalization is 
applied to the system up to 2N = 20 sites. Checking numerically that magnetic states with 
the magnetization per site m obey a conformal field theory with conformal anomaly c = 1 
for 1/4 < m < 3/4, we use the finite-size scaling of the conformal invariance to obtain a 
magentization curve in the thermodynamic limit. In the magnetizatin curve a plateau appears 
at m = 1/4. We also calculate two critical exponents r\ and rf for 1/4 < m < 3/4, which control 
the asymptotic behavior of the transverse and parallel spin correlation functions. We check the 
relation r\rf = 1, which universally holds for a c = 1 conformal field theory. 
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Recently, spin chains with an alternating array of dif- 
ferent spins have attracted some current interest. The- 
oretically, Lieb and Mattis showed that the alternat- 
ing spin antiferromagnetic Heisenberg chains with some 
next-nearest-neighour interactions have a ground state 
with large total spin. ]l[] Very recently, numerical calcu- 
lations of numerical diagonalization, [^|,^] density-matrix 
renomalization method Q and quantum Monte Carlo 
method § have been applied to the alternating spin 
Heisenberg chains. Experimentally, the alternating spin 
chains have been found as quasi-one-dimensional ferri- 
magnetic chains. || On the other hand, the integrable al- 
ternating spin models have been constructed and solved 
via Bethe ansatz. M Although its Hamiltonian has 
complicated nearest- and next-nearest-neighbour inter- 
actions, the integrable model has a single ground state. 

Generalizing the Lieb-Shultz-Mattis theorem, || Os- 
hikawa et al gave a presence condition of a magnetic 
plateau for spin chains in a magnetic field, j^] The mag- 
netic plateau has been observed in some spin models; 
for examples, an S = 1/2 antiferromagnetic chain with 
period 3 exchange coupling |i"o| ] and an S = 1 antiferro- 
magnetic chain with alternating bond. |nj The alternat- 
ing spin chain is one of simple examples for the magnetic 
plateau. Sakai and Takahashi performed the numercal 
calculation for S = 1 antiferromagnetic Heisenberg chain 
in a magnetic field and revealed its magnetic and criti- 
cal properties in the thermodynamic limit by using the 
technique of the conformal field theory. jl2| Applying 
the Bethe ansatz to the integrable altenating spin chain 
with 5=1/2 and 1, Fujii et al calculated a magnetiza- 
tion curve and some critical exponents. |D| They did not 
observe any magnetic plateau but a strange cusp in the 
magnetization curve. In this paper, following the pro- 
cedure developed in ref.12, we study the magnetic and 
critical properties of the alternating spin antiferromag- 



netic Heisenberg chain with S — 1/2 and 1 in a magnetic 
field at T = by numerical diagnoalization and use the 
finite-size scaling of the confomal field theory to find out 
its magnetic and critical properties in the themodynamic 
limit. 

The alternating spin chain of 2N sites consists of N 
spins cri,cr3,- ■ •, (T2N-1 of spin 1/2 and N spins S^2, SU, 
■ ■ •) &2N of spin 1. This spin chain in a magnetic field H 
is described by the Hamiltonian 



H, = Tip, 
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Ho — 2i-l ' S'2i + S2i ■ T2i+l) 
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where we imopse periodic boundary conditions; (T2N+1 = 
<Ti and S2N+2 = ^2- The Hamiltonian TL is invariant un- 
der two-site translation and rotation about z axis. Thus, 
all eigenstates of the Hamiltonian Ti can be classified 
by the magnetization M — J2iLi ( a 2i-i + ^ii) an d the 
wave vector k (k = 2im/N, n = 0, 1, 2, • • • , N - 1 ). We 
calculate the lowest energy Ek (2 AT, M) of the Hamilto- 
nian 7in in the subspace where the system has the mag- 
netization M and the wave vector k, by the numerical 
diagonalization with Lanczos algorithm up to 2N = 20. 
We define E(2N, M) as the lowest energy of the Hamil- 
tonian 7i in the subspace with the magnetization M. In 
our case, we find that E(2N, M) = E k=0 (2N, M). 

The conformal field theory predicts that if the lowest- 
energy state with the magnetization per site m is mass- 
less, the size-dependence of the lowest energy per site has 
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the form juj 

(27V -> oo) 

where we must vary 27V and M with the magnetization 
per site fiexd at m — M/(2N). e(m) is the lowest en- 
ergy per site and v s (m) is the sound velocity. They de- 
pend on the magnetization per site m. c is the conformal 
anomaly. 



Thus, the sound velocity v s (m) is estimated by 

v s {m) = ^ (E 2 „ /N (2N, M) - E(2N, M)) (6) 

for the largest 2N . The results of the estimation for the 
conformal anomaly c is shown in Table I. The value of 
c for m — 7/20, 3/8 and 5/12 is slightly larger than 1. 
As we will see later, the universal relation between two 
exponents r\ and rf for a c = 1 conformal field theory 
holds in our case. Thus we conclude that c = 1 for 
1/4 < m < 3/4. 

We consider a magnetization curve of the alernating 
spin chain in two regions of < m < 1/4 and 1/4 < 
m < 3/4. For < H < H c ±, where the value H c \ is 
equal to the energy gap A, the lowest-energy state of 
the Hamiltonian 7i has the magnetization M = N/2. 
This means that m = 1/4 for < H < H cX {= A). For 
1/4 < to < 3/4, owing to the confomal invariance, the 
spin-excitation energy of the Hamiltonian TL depends on 
the system size 2N as 

E{2N, M + 1) - E{2N, M) - H ~ nv s r]^- (7) 



Fig. 1. Plots of E(2N,M)/(2N) vs 1/(27V) 2 with m = M/(2N) 
fixed at to = 1/4, 5/12, 1/2 and 7/12. The origin is shifted 
along the vertical axis without changing the scale. The values 
of points A and B are as follows; [A, B] = [-0.721,-0.771], 
[-0.400, -0.450], [-0.195, -0.245] and [-0.040, -0.090] for to = 
1/4, 5/12, 1/2 and 7/12, respectively. The solid lines show the 
fitting of two values E(2N, M)/(2N) for the largest and next- 
largest 2N to the form (5). The dotted line is guide to the eye. 



E(2N, M) - E(2N, M-l)-H~ ~™ s r,^ (8) 

in the limit 2N — > oo, [Q where we vary 2N and 
M with m = M/(2N) fixed. 77 is a critical exponent 
of correlation function for spin-excitations, which we 
will estimate in the following. In Fig. 2 , we plot 



In Fig.l, we show plots of E(2N,M)/(2N) versus 
l/(2iV) 2 for to = 1/4, 5/12, 1/2, 7/12. The plots are 
alomst linear for 777 = 5/12, 1/2, 7/12 but values for 
777 = 1/4 converge fast than 1/(2N) 2 . Since the ground 
state of the alternating spin chain has the total spin N/2 
and there is an energy gap A between the ground state 
and the lowest-energy state of the total spin N/2 + 1,0] 
they suggets that the lowest energy state is massless for 
1/4 < m < 3/4 but massive for < m < 1/4. ( We use 
the energy gap A = 1.767 ± 0.003 obtained in ref.5.) In 
the following, assuming that the lowest-energy state is 
massless for 1/4 < 777 < 3/4, we apply the scaling law 
of conformal invariance to find out magnetic and criti- 
cal properties in the thermodynamic limit. To determine 
the conformal anomaly c, we need two values, the sound 
velocity v s and the gradient of the plot E(2N, M)/{2N) 
versus l/(2A r ) 2 , — ircv s (m)/6. To estimate the gradient, 
we use two values of E(2N, M) for the largest and next- 
largest 27V such that the magnetization per site is fixed 
at 777 = M/(2N), up to 27V = 20. The next-lowest- 
energy state appears at the wave vector k = 2-k/N . 
The gap between the lowest- and next-lowest energies, 
E 2w /n{2N,M) - E(2N,M), has the dependence on the 
system size as 

E 2 „ /N (2N, M) - E(2N, M) ~ ^ s (m) (5) 

(27V -► 00). 



Fig. 2. Plots versus 1/(27V) of spin-excitation gaps, E(2N, M + 
1) - E(2N, M) (open circles) and E(2N, M) - E(2N, M - 1) 
(crosses), with m = M/(2N). The solid lines show the fitting of 
two values of E(2N, M+1)~E{2N, M) ( E(2N, M) — E(2N, M — 
1) ) for the largest and next-largest 2N to the form (7) ( (8) ). 
For to = 1/4, the gap E(2N, M + 1) - E(2N, M), connected by 
a dotted line, converges to a finite value A. 

E{2N,M+1)-E(2N,M) &ndE(2N,M)-E(2N,M-l) 
versus 1/(27V). The plots are almost linear at least 
for 777 = 5/12, 1/2 and 7/12. For m = 1/4, the gap 
£;(27V, M+1)-E(2N, M) converges faster than 1/(27V). 
Combinning eqs. (5) and (7) (or (8)) and taking the 
limit 27V — > 00 , we get the relation between the magen- 
tic field H and magnetization per site 777, e'(m) = H, 
which gives us a magnetization curve at T — 0. In the 
following, we estimate e'(m) for several m's and con- 
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m 7/20 3/8 5/12 9/20 1/2 11/20 7/12 5/8 13/20 



c 1.20 1.23 1.18 1.07 1.07 1.04 1.08 1.07 1.04 



Table I. Conformal anomaly c estimated from egs. (5) and (6). 



nect them by a smmoth curve. To etsimate e'(m) for 
1/4 < to < 3/4, we use the largest- and next-largest-size 
data for E(2N, M + 1) - E(2N, M) in the form (7) and 
do the same treatment for E(2N, M) ~ E(2N, M - 1) in 
the form (8). Two results of e'(m) for 9/20 < to coin- 
cide with each other within a difference less than 0.5%. 
The difference for m = 7/20, 3/8 and 7/12 is less than 
1%. Thus we regard the averge of the two results as the 
extraporated value of e'(m). Using the extrapolated val- 
ues of e'(m), we draw in Fig. 3 the magnetization curve at 
T = in the thermodynamic limit. We see that a plateau 



whose definition will be given later. First we consider 
the critical exponent r\. The exponent r\ appears in the 
asymptotic forms (7) and (8) of the spin-excitation en- 
ergy gaps. Thus, substitutting eq.(8) from eq.(7) and 
using eq.(6), we get the relation 

V = 

E(2N, M + 1) + E(2N, M — 1) — 2E(2N, M) 
E 2v/N (2N,M)-E(2N,M) 

(9) 

The results of t] evaluated by eq. (9) for 2N = 16, 18 
and 20 are plotted in Fig. 4 versus magnetization per site 
m. We determine r\ in the limits of m — > (1/4)+ and 



Fig. 3. Magnetization curve at the thermodynamic limit. Open 
circules show plots of magnetization per site m versus H(= 
e'(m)) which is estimated by averaging two results extrapolated 
by eqs. (7) and (8). A plateau appears at m = 1/4. 



Fig. 4. Critical exponent r) estimated from eq. (9) for 27V = 16, 
18 and 20. 



appears at m = 1/4. The presence of this plateau is ex- 
pected from the Lieb-Mattis theorem, Q which shows 
that the ground state of the alternating spin antiferro- 
magnetic Heisenberg chain has the total spin N/2 and 
has an energy gap between the ground state and the 
lowest-energy state of spin N/2 +1. It is also consis- 
tent with the generalized LSM theorem proposed by Os- 
hikawa et al Q. The generalized LSM theorem can be 
applied to general spin chains with axial symmetry. The 
quantization condition of a magnetic plateau in our case 
is that 1/2 + 1 — 2m =intger, which holds for to = 1/4. 

Now we calculate two critical exponents, r\ and rf : , for 
correlation functions of the transverse and parallel com- 
ponents to a magnetic field, which behave as (— l) r r~ n 
for the transverse component and r~'' cos(2fci?r) for the 
parallel one, where r is the distance between two opera- 
tors in the correlation function and 2k F is a wave vector 



m — » (3/4) — , by extrapolating the value of 77 for M = 
N/2 + 1 (to = (1/4)+) and M = 3N/2-1 (to = (3/4)-) 
linearly to l/(2iV) for 2N = 18 and 20. The results are 
77 = 0.460 for m = 1/4 and 77 = 0.493 for m = 3/4. 

Next we consider the critical exponent ?7 Z for the spin- 
spin correlation function. In the numerical calculation of 
E k (2N,M) for k = 0, ir/N, 2n/N,- ■ 2tt(1 - l/N), we 
find that a soft mode exits at k — 2txM /N — tt = 2kp for 
M < N and k = 3?r- 2nM/N = 2k F for M < N . If the 
soft mode become a gapless excitation as 2N —> 00, the 
conformal invariance predicts that its energy gap from 
the lowest energy depends on the system size 27V as 

E 2kF (2iV, M) ~ E{2N, M) ~ nv^ ^ (10) 

in the limit 2N — » 00. |L5| We numerically check this 
asymptotical behavior of the energy gap for to = 5/12, 



1 



1/2 and 7/12. Thus we assume that the gapless excita- 
tion exits at k = 2kp in the infinite system. To estimate 
r] z , we use the relation 

z _ E 2kF (2N,M)-E(2N,M) 
1 E 2n/N (2N,M)-E(2N,M) [ ' 

derived from eqs. (6) and (10). The results of rf for 
27V = 16, 18 and 20 are plotted in Fig.5. 
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0.617 


2.000 


1.234 
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0.608 


1.825 


1.109 
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0.606 


1.744 


1.056 
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0.601 


1.721 


1.034 
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0.592 


1.717 


1.016 
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0.578 


1.777 


1.027 


12 


0.563 


1.833 


1.031 


13 


0.548 


1.902 


1.043 


14 


0.539 


2.000 


1.078 



Table II. Exponents r\ and rf estimated from eqs. (9) and (11) 
for 2N = 20 and the values of r\r\ z . 



Fig. 5. Critical exponent rj z estimated from eq. (11) for 2N = 16, 
18 and 20. 

We compare the above results of rj and rf with those 
for 5 — 1/2 and S = 1 antiferromagnetic Heisenberg 
chains. Let us first note the behaviors of rj and r] z for 
5 = 1/2 and 5 = 1. For 5 = 1/2, when m varies from 
to 1/2, rj monotonously decreases from 1 to 1/2 and 
if monotonously increases from 1 to 2. |l6| For 5=1, 
rj = 1/2 and rf — 2 at to — and m — 1. rj and rf have 
a minimum and a maximum near to — 1/3, respectively. 
]l2] | The parallel spin correlation for 5 = 1/2 is stronger 
than that for 5=1. In the alternating spin chain, 77 and 
rf have a maximum and a minimum, respectively. The 
parallel spin correlation for the alternating spin chain is 
stronger than that for 5=1 but the transverse one is 
weaker. The alternating spin is similar in this aspect to 
5 = 1/2. 

Finally, let us check that two estimated critical expo- 
nents 77 and rf satisfy an univeral relation r\rf = 1 for 
1/4 < to < 3/4. This relation holds universally for a 
c = 1 conformal field theory. |i~7| The values of r\rf for 
each 77i in the case of 27V = 20 are shown in Table II. It 
shows that the relation r\rf = 1 is satisfied within errors. 
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